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ABSTRACT 

The BL Lacertae object 3C 66A was detected in a flaring state by the Fermi Large 
Area Telescope (LAT) and VERITAS in 2008 October. In addition to these gamma-ray 
observations, F-GAMMA, GASP-WEBT, PAIRITEL, MDM, ATOM, Swift, and Chan- 
dra provided radio to X-ray coverage. The available light curves show variability and, 
in particular, correlated flares are observed in the optical and Fermi-L AT gamma-ray 
band. The resulting spectral energy distribution can be well fit using standard leptonic 
models with and without an external radiation field for inverse-Compton scattering. It 
is found, however, that only the model with an external radiation field can accommodate 
the intra-night variability observed at optical wavelengths. 

Subject headings: BL Lacertae objects: individual (3C 66A) — galaxies: active - 
gamma rays: observations 



Introduction 



The radio source 3C 66 ([Bennett! Il962j ) was shown by iMackavl (|197ll ) and iNorthoverl (|1973l ) 
to actually consist of two unrelated radio sources separated by 0.11° : a compact source (3C 66A) 
and a resolved galax y (3C 66B). 3C 66A was subsequen tly identified as a q uasi-stellar object by 
Wills k, Wills! (|1974j ). and as a BL Lacertae object by iSmith et al.l (|1976l ) based on its optical 
spectrum. 3C 66A is now a well-known blazar which, like other active galactic nuclei (AGN), is 
thought to be powered by accretion of material ont o a supermassive black hole located in the central 
region of the host galaxy (jUrry &: Padovanilll995l ). Some AGN present strong relativistic outflows 
in the form of jets, where particles are believed to be accelerated to ultra-relativistic energies and 
gamma rays are subsequently produced. Blazars are the particular subset of AGN with jets aligned 
to the observer's line of sight. Indeed, the jet of 3C 66A has been imaged using very long baseline 



interferometry (VLBI) ([Taylor et al.lll996l ; iJorstad et al 



20071 ) and superluminal motion has been inferred (jJorstad et al 



2001 



Marscher et al 



2001 



2002 



Britzen et al. 



Britzen et al 



2008). This is 



indicative of the relativistic Lorentz factor of the jet and its small angle with respect to the line of 
sight. 

BL Lacs are known for having very weak (if any) detectable emission lines, which makes 
de termination of thei r redshift quite difficult. The redsh ift of 3C 66A was r eported as z = 0.444 
bv lMiller et al.l (|1978l ) and also (although tentatively) bv lKinnev et al.l (|199ll ). Each measurement 



120 Now at DESY, Platanenallee 6, 15738 Zeuthen, Germany 

121 Now at Institut fair Physik und Astronomie, Universitat Potsdam, 14476 Potsdam-Golm, Germany; DESY, Pla- 
tanenallee 6, 15738 Zeuthen, Germany 

122 Now at Los Alamos National Laboratory, MS H803, Los Alamos, NM 87545 
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however, is based on the measurement of a single line and is not reliable (jBramel et al.ll2005l ). 
Recent efforts (described in Section 12, 5p to provide further constraints have proven unsuccessful. 

Similar to other blazars, the spectral energy distribution (SED) of 3C 66A has two pro- 
nounced peaks, which sugges ts that at least two different physical emission processes are at work 
(je.g.. Joshi &: Bottcherl 120071 ). The first peak, extending from radio to soft X-ray frequencies, is 
likely due to synchrotron emission from high-energy electrons, while different emission models have 
been proposed to explain the second peak, which extends up to gamma-ray energies. Given the 
location of its synchrotron peak ( < 1Q 15 Hz) ; 3C 66 A is further sub-classified as an intermediate 
synchrotron peaked blazar (ISP) (jAbdo et al.ll2010d ). 



The models that have been proposed to explain gamma-ray emission in blazars can be roughly 
categorized into leptonic or hadronic, depending on whether the accelerated particles responsible 
for the gamma-ray emission are primarily electrons and positrons (hereafter "electrons") or protons. 
In leptonic models, high-energy electrons produce gamma rays via inverse Compton scattering of 
low-energy photons. In synchrotron self-Compton (SSC) models, the same population of electrons 
responsible for the observed gamma rays generates the low-energy photon field through synchrotron 
emission. In external Compton (EC) models the low-energy photons originate outside the emission 
volume of the gamma rays. Possible sources of target ph otons include: accretion-disk photons 
radiated directly into the jet (jDermer &: Sch lickeiscr 1993), accretion-disk photons scattered by 



emission-line clouds or dust into the jet (ISikora et al 



back into the jet by broad-line emission cl ouds (iGhisellini &: Madaulll996l ). jet emission from an 



1994 ), synchrotron radiation re-scattered 



outer slow jet sheet (IGhisellini et al.ll2005l ). or emission from faster or slower portions of the jet 
(jGeorganopoulos fc Kazanasll2004l ). In hadronic m odels, gamma rays are produced by high-energy 



protons, either via proton synchrotron radiation (IMiicke et al.l 120031). or via secondary emission 



Bottcher 



(|2007l ) and references therein 



from photo-pion and photo-pair-production reactions (see 
for a review of blazar gamma-ray emission processes) . 

One of the main obstacles in the broadband study of gamma-ray blazars is the lack of si- 
multaneity, or at least contemporaneousness, of the data at the various wavelengths. At high 
energies the situation is made even more difficult due to the lack of objects that can be detected by 
MeV/GeV and TeV observatories on comparable time scales. Indeed, until recently the knowledge 
of blazars at gamma-ray energies had been obtained from observations performed in two disjoint 
energy regi mes: i) the high energy (HE) range (20 MeV< E < 10 GeV), studied in the 1990s 
by EGRET (jThompson et al.lll993h . a nd ii) the very -high-energy (VHE) regime (E > 100 GeV) 
observed by ground-based instruments ( Weekes 2000 ). Onl}Q Mark arian 421 was detect ed by both 
EGRET and the first imaging atmospheric Cherenkov telescopes (jKerrick et al.lll995l ). Further- 
more, blazars detected by EGRET at MeV/GeV energies are predominantly flat-spectrum radio 
quasars (FSRQs), while TeV blazars are, to date, predominantly BL Lacs. It is important to un- 
derstand these observational differences since they are likely related to the physics of the AGN 



1 Markarian 501 was marginally detected by EGRET only during a few months in 1996 1 Kataoka et al. 



1999l l. 
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j Cavaliere fc D'Elial 120021 ) , or to the evolution of blazars over cosmic time ([Bottcher Dermer 
20021 ) . 



The current generation of gamma-ray instruments (AGILE, Fermi, H.E.S.S., MAGIC and 
VERITAS) is closing the gap between the two energy regimes due to improved instrument sensi- 
tivities, leading us towards a deeper and more c omplete characterization of blazars as high-energy 
sources and as a population (jAbdo et al.ll2009bl ). An example of the successful synergy of space- 
borne and ground-based observatories is provided by the joint observations of 3C 66A by Fermi and 
VERITAS during its strong flare of 2008 October. The flare was originally reported by VERITAS 
(|Swordvll2008l ; lAcciari et al.l I2009P ) , and soon after, con temporaneous variability was also det ected 
at op tical to infrared wavelengths (jLarionov et al.ll2008l ). and in the Fermi-L AT energy band (jTosti 
2008). Follow-up observations were obtained at radio, optical and X-ray wavelengths in order to 
measure the flux and spectral variability of the source across the electromagnetic spectrum and to 
obtain a quasi-simultaneous SED. This paper reports the results of this campaign, including the 
broadband spectrum and a model interpretation of this constraining SED. 



2. Observations and Data Analysis 
2.1. VERITAS 



The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an array of 
four 12m diameter im aging atmospheric Cherenkov telescopes (IACTs) in southern Arizona, U.S.A. 
(|Acciari et al.1 l2008bh . 3C 66A was observed with VERITAS for 14 hours from 2007 September 
through 2008 January and for 46 hours between 2008 September and 2008 November. These 
observations (hereafter 2007-2008 data) add up to ~32.8 hours of li ve time after data q uality 



selection. The data were analyzed following the procedure described in lAcciari et al.l (|2008bl ). 



As reported in lAcciari et al.l (|2009l ). the average spectrum measured by VERITAS is very soft, 



yielding a photon index V of 4.1 ±0A s t a t ± 0.6 S ^ S when fitted to a power law dN/dE oc E 
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cm s , 



The average integral flux above 200 GeV measured by VERITAS is (1.3 ± 0.1) x 10" 
which corresponds to 6% of the Crab Nebula's flux above this threshold. In addition, a strong 
flare with night-by-night VHE-flux variability was detected in October 2008. For this analysis the 
VERITAS spectrum is calculated for the short time interval October 8 - 10, 2008 (MJD 54747- 
54749; hereafter flare period), and for a longer period corresponding to the dark rurQ where most 
of the VHE emission from 3C 66A was detected (MJD 54734 - 54749). It should be noted that the 
flare and dark run intervals overlap and are therefore not independent. Table [U lists the relevant 
information from each data set. 



2 Imaging atmospheric Cherenkov telescopes (IACTs) like VERITAS do not operate on nights with bright moon- 
light. The series of nights between consecutive bright-moonlight periods is usually referred as a dark run. 
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As shown in Figure [TJ the flare and dark run spectra are very soft, yielding nearly identical 
photon indices of 4.1 ±0.6 sta t ±0.6 S2/S , entirely consistent with that derived from the full 2007-2008 



-11 —2 

cm s 



data set. The integral flux above 200 GeV for the flare period is (2.5 ± 0.4) x 10" 
while the average flux for the dark run period is (1.4 ± 0.2) x 10 -11 cm -2 s _1 . The extragalactic 
background light (EBL) de-absorbed spe ctral points for the dark run calculated using the optical 
depth values of (jFranceschini et al.ll2008l ) and assuming a nominal redshift of z = 0.444, are also 
shown in Figured! These points are well fit by a power law function with T = 1.9 ± 0.5. 



2.2. Fermi-L AT 



The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope is a pair- 
conversion detector s ensitive to gamma rays with energies between 20 MeV and several hundred GeV 
(jAtwood et al.ll2009i ). Since launch the instrument has operated almost exclusively in sky survey 
mode, covering the whole sky every 3 hours. The overall coverage of the sky is fairly uniform, 
with exposure variations of < 15% around the mean value. The LAT data are analyzed using 
ScienceTools v9rl5p5 and instrument response functions P6V3 (available via the Fermi science 
support centejfl). Only photons in the diffuse event class are selected for this analysis because of 
their reduced charged-particle background contamination and very good angular reconstruction. A 
zenith angle <105° cut in instrument coordinates is used to avoid gamma rays from the Earth limb. 
The diffuse emission from the Galaxy is modeled using a spatial model (gll_iem_v02.f it) which 
was refined with Fermi-LAT data taken during the first year of operation. The extragalactic diffuse 
and residual instrumental backgrounds are modeled as an isotropic component an d are included 
the fhfl The data are analyzed with an unbinned maximum likelihood technique (jMattox et al 



m 



19961 ) using the likelihood analysis software developed by the LAT team. 



Although 3C 66A was detected by EGRET as source 3E G J0222+4253 (|Hartman et alJll999h . 
detailed spatial and timing analyses by iKuiper et al.l (|2000l ) showed that this EGRET source ac- 
tually consists of the superposition of 3C 66A and the nearby millisecond pulsar PSR J0218+4232 
which is 0.96° distant from the blazar. This interpretation of the EGRET data is verified by 
Fermi-L AT , whose improved angular resolution permits the clear separation of the two sources 
as shown in Figure [2] Furthermore, the known pulsar period is detected with high confidence in 
the Fermi-LAT data (jAbdo et al.l l2009al ) . More importantly for this analysis, the clear separation 
between the pulsar and the blazar enables studies of each source independently in the maximum 
likelihood analysis, and thus permits an accurate determination of the spectrum and localization 
of each source, with negligible contamination. 



Figure [2] also shows the localization of the Fermi and VERITAS sources with respect to blazar 



3 http://fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/overview. html 
4 http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 
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3C 66A and radio galaxy 3C 66B (see caption in Figure [2] for details). It is clear from the map 
that the Fermi-LAT and VERITAS localizations are consistent and that the gamma-ray emission 
is confidently associated with the blazar and not with the radio galaxy. Some small contribu- 



tion in the Fermi-LAT data from radio galaxy 3C 66B as suggested by lAliu et al.l (|2009l ) and 
Tavecchio fc Ghisellinil (|2009h cannot be excluded given the large spill-over of low-energy photons 
from 3C 66A at the location of 3C 66B. This i s due to the l o ng ta ils of the Fermi-LAT point- 
spread- function at low energies as described in lAtwood et al.l (120091 ) . Nevertheless, considering 
only photons with energy E > 1 GeV, the upper limit (95% confidence level) for a source at the 
location of 3C 66B is 2.9 x 10~ 8 cm — 2 s -1 for the dark run period (with a t est statistic^ TS = 1.3). 
For the 11 months of data corresponding to the first Fermi-LAT catalog (jAbdo et al.ll2010al ) the 
upper limit is 4.9 x lO^cm^s -1 (TS = 5.8). 

As in the analysis of the VERITAS observations, the Fermi-LAT spectrum is calculated for the 
flare and for the dark run periods. The Fermi flare period flux F(E >100MeV) = (5.0 ± 1.4 stat ± 
0.3 sys ) x 10 cm - 2 s is consistent within errors with the dark run flux of (3.9 ± 0.5 s t a t =t 0.3 sys ) 
x 10" 7 



cm 



. In both cases the Fermi-LAT spectrum is quite hard and can be described by a 
power law with a photon index T of 1.8 ± 0.1 s t a t ± 0.1 sys and 1.9 ± 0.1 stat ± 0.1 sys in the flare period 
and dark run intervals, respectively. Both spectra are shown in the high-energy SED in Figure [TJ 



2.3. Chandra 



3C 66A was observed by the Chandra observatory on October 6, 2008 for a total of 37.6 ksec 
with the Advanced CCD Imaging Spectrometer (ACIS), covering the energy band between 0.3 and 
10 keV. The source was observed in the continuous clocking (CC) mode to avoid pile-up effects. 
Standard analysis tools (CIAO 4.1) and calibration files (CALDB v3.5.0) provided by the Chandra 
X-ray centejfl are used. 

The time-averaged spectrum is obtained and re-binned to ensure that each spectral channel 
contains at least 25 background-subtracted counts. This co ndition allows the use of the x 2 quality- 
of-fit estimator to find the best fit model. XSPEC vl2.4 (lArnaudlll996l ) is used for the spectral 
analysis and fitting procedure. 



Two spectral models have been used to fit the data: single power law and broken power 
law. Each model includes galactic Hi column density ( N#Q a i = 8.99 x 10 20 cm~ 2 ) according to 
Dickey fc Lockman (1990), where the photoelectric absorption is set with the XSPEC model p/ia&J 7 ]. 



5 The test statistic value (TS) quantifies the probability of having a point source at the location specified. It is 
roughly the square of the significance value: a TS of 25 corresponds to a signal of approximately 5 standard deviations 



(Abdo et al. 2010a) 



6 http://cxc. harvard.edu/ciao/ 

7 http: / /heasarc. gsfc.nasa.gov / docs/software/lheasoft / xanadu / xspec / manual/XSmodelPhabs.html 
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Interval 


Live Time [hr] 


iVon 


Nob 


Alpha 


Excess 


Significance [a] 


flare 


6.0 


1531 


7072 


0.121 


678.3 


18.0 


dark run 


21.2 


3888 


20452 


0.125 


1331.5 


22.2 


2007-2008 


28.1 


7257 


31201 


0.175 


1791 


21.1 



Table 1: Results from VERITAS observations of 3C 66A. Live time corresponds to the effective 
exposure time after accounting for data quality selection. A r o n (-ZVoff) corresponds to the number 
of on (off )-source events passing background-rejection cuts. Alpha is the normalization of off-source 
events and the excess is equal to Aon — oiNob ■ The signi ficance is e xpres sed in number of st a ndard 
deviations and is calculated according to equation (17) of lLi fc Mai (|1983l ). See lAcciari et al.l (j.2009) 
for a complete description of the VERITAS analysis. 
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Fig. 1 . — Gamma-ray SED of 3C 66A including Fermi-LAT and VERITAS data for the flare (red symbols) and 
dark run (blue symbols) intervals. The Fermi-LAT spectra is also shown here as "butterfly" contours (solid lines) 
describing the statistical error on the spectrum ( Abdo et al. 2009b). The previously reported Fermi-LAT 6-month- 



average spectrum IjAbdo et al.ll2010bf ) is also shown here (green circles) and is lower than the spectrum obtained 
during the campaign. The average 2007-2008 VERITAS spectrum originally reported in Acciari et al. 1 20091 ) is 
displayed with green triangles. In all cases the upper limits are calculated at 95% confidence level. The de-absorbed 
dark run spectra obtained using the optical depth values of Franceschini et al. ( 2008h are also shown as open circles 
and open squares for redshifts of 0.444 and 0.3, respectively. 
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Right Ascension [h] 




Counts per pixel 



Fig. 2. — Smoothed count map of the 3C 66 A region as seen by Fermi-L AT between September 
1, 2008 and December 31, 2008 with E>100 MeV. The color bar has units of counts per pixel and 
the pixel dimensions are 0.1° x 0.1°. The contour levels have been smoothed and correspond to 
2.8, 5.2, and 7.6 counts per pixel. The locations of 3C 66A and 3C 66B (a radio galaxy that is 
0.11° away) are shown as a cross and as a diamond, respectively. The location of millisecond pulsar 
PSR 0218+4232 is also indicated with a white cross. The magenta circle represents the VERITAS 
localization of the VHE source (RA; DEC) = (2 h 22 m 41.6 s 



±1.7 s sta t ±6.0 s sys ; 43° 02' 35.5" 



± 21" s tat i 1'30" sys ) as reported in lAcciari et al.l (|2009l ). The blue interior circle represents the 
95% error radius of the Fermi-LAT localization (RA; DEC) = (02 h 22 m 40.3 s ± 4.5 s ; 43° 02' 18.6" 



± 42.1") as reported in the Fermi-LAT first source catalog ( Abdo et al.ll2010al ). All positions are 
based on the J2000 epoch. 
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An additional local H I column density was also tried but in both cases the spectra was consistent 
with pure galactic density. Consequently, the column density has been fixed to the galactic value 
in each model and the obtained results are presented in Table [2J An F-test was performed to 
demonstrate that the spectral fit improves significantly when using the extra degrees of freedom of 
the broken power law model. Table [2] also contains the results of the F-test. 



2.4. Swift XRT and UVOT 



Following the VERITAS detection of VHE emission from 3C 66A, Target of Opportunity 
(ToO) observations of 3C 66A with Swift were obtained for a total duration of ~10 ksec. The Swift 
satellite observatory co mprises an UV-Opti cal telescope (UVOT), an X-ray telescope (XRT) and a 
Burst Alert Telescope (jGehrels et al.ll2004l ). Data reduction and calibration of the XRT data are 
performed with HEASoft v6.5 standard tools. All XRT data presented here are taken in Photon 
Counting (PC) mode with negligible pile-up effects. The X-ray spectru m of each observation is fi t 
with an absorbed power law using a fixed Galactic column density from lDickev fc Lockmanl ([1990), 
which gives good x 2 values for all observations. The measured photon spectral index ranges between 
2.5 and 2.9 with a typical statistical uncertainty of 0.1. 

UVOT obtained data through each of six color filters, V, B and U together with filters defining 
three ultraviolet passbands, UVW1, UVM2 and UVW2 with central wavelengt hs of 260 nm, 220 
nm and 193 nm respectively. The data are calibrated using standard techniques (IPoqle et al.ll2 008) 



Schlegel et al 



and c orrected for Galactic extinction by interpolating the absorption val ues from 
(j 19981 ) (Eb-v = 0.083 mag) with the galactic spectral extinction model of iFitzpatrickl (|1999i ). 



2.5. Optical to Infrared Observations 



The R magnitude of the host galaxy of 3C 66 A is ~19 in the optical band (jWurtz et al 



19961 ). Its contribution is negligible compared to the typical AGN magnitude of R <15, therefore 



host-galaxy correction is not necessary. 

GASP-WEBT: 3C 66A is continuously monitored by telescopes affiliated to the GLAST 



AGILE supp ort program of the Whole Earth Blazar Telescope (GASP-WEBT; see IVillata et al 



(|2008l . l2009i )). These observations provide a long-term light curve of this object with complete 
sampling as shown in Figure El During the time interval in consideration (M JD 54700 - 54840) , 
several observatories (Abastumani, Crimean, L'Ampolla, Lulin, New Mexico Skies, Roque de los 
Muchachos (KVA), Rozhen, Sabadell, San Pedro Martir, St. Petersburg, Talmassons, Teide (BRT), 
and Tuorla) contributed photometric observations in the R band. Data in the J, H and K band 
were taken at the Campo Imperatore observatory. A list of the observatories and their location is 
available in Table [3l 
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Single Power Law Model 


r 


Flux [10 12 ergs cm 2 s 1 ] 


X 2 /d.o./. 


2.99±0.03 


3.47±0.06 


1.21 (232.6/193) 


Broken Power Law Model 




r 2 


Flux [10 12 ergs cm 2 s '] 


Break [keV] 


X 2 /d.o.f. 


F-tcst Prob. 


3.oa±S;S 


2 24+ - 23 


O ro+0.07 

°- d °-0.08 


o n+0.5 
°-°-0.3 


0.97 (185.2/191) 


3.47 x 1(T 10 



Table 2: Best-fit model parameters for a fit performed to the Chandra data in the 1 — 7 keV energy 



range. The galactic Nffr?n./ value is fixed to 8.9 9 x 10 20 cm z , the value of the galactic Hi column 



density according to lDickev &: Lockmanl (|199Cll ). Errors indicate the 90% confidence level. 



Observatory 


Location 


Web page 




Radio Observatories 


Grimcem r\,a,aio Obs. 


Ukraine 


www.crao.crimea.ua 


Effelsberg 


Germany 


www.mpifr.de / english / radiotelescope 


TD A Al 

IKAlvl 


Spain 


www.iram-institute.org/EN/30-meter-telescope.php 


ivieuiciiia 


Italy 


www.med.ira.inaf.it 


Metsahovi 


Finland 


www.metsahovi.fi/en 


Noto 


Italy 


www.noto.ira.inaf.it 


UMRAO 


Michigan, USA 


www.astro.lsa.umich.edu/obs/radiotel 




Infrared Observatories 


Campo Imperatore 


Italy 


www.oa-teramo.inaf.it 


PAIRITEL 


Arizona, USA 


www.pairitel.org 




Optical Observatories 


Abastumani 


Georgia 


www . genao . org 


Armenzano 


Italy 


www.webalice.it /dcarosati 


ATOM 


Namibia 


http: / /www.lsw.uni-heidelberg.de /projects/hess/ ATOM/ 


Crimean Astr. Obs. 


Ukraine 


www.crao.crimea.ua 


Kitt Peak (MDM) 


Arizona, USA 


www.astro.lsa.umich.edu/obs/mdm 


L'Ampolla 


Spain 




Lulin 


Taiwan 


www . lulin . ncu . edu . tw/ english 


New Mexico Skies Obs. 


New Mexico, USA 


www.nmskies.com 


Roque (KVA) 


Canary Islands, Spain 


www .otri.iac.es/eno/nt.htm 


Rozhen 


Bulgaria 


www.astro.bas.bg/rozhen.html 


Sabadell 


Spain 


www.astrosabadell.org/html/es/observatoriosab.htm 


San Pedro Martir 


Mexico 


www.astrossp.unam.mx/indexspm.html 


St. Petersburg 


Russia 


www.gao.spb.ru 


Talmassons 


Italy 


www.castfvg.it 


Teide (BRT) 


Canary Islands, Spain 


www.telescope.org 


Torino 


Italy 


www.to.astro.it 


Tuorla 


Finland 


www.astro.utu.fi 


Valle d' Aosta 


Italy 


www.oavda.it / english / osservatorio 




Gamma- 


ray Observatories 


VERITAS 


Arizona, USA 


veritas.sao.arizona.edu 



Table 3: List of ground-based observatories that participated in this campaign. 
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MDM: Following the discovery of VHE emission, 3C 66A was observed with the 1.3m telescope 
of the MDM Observatory during the nights of Oct. 6 - 10, 2008. A total of 290 science frames in 
U, B, V, and R bands (58 each) were taken throughout the entire visibility period (approx. 4:30 - 
10:00 UT) during each night. The light curves, which cover the time around the flare, are presented 
in Figured! 

ATOM: Optical observations for this campaign in the R band were also obtained with the 0.8 
m optical telescope ATOM in Namibia which monitors this source periodically. Twenty photometric 
observations are available starting on MJD 54740 and are shown in Figures [3] and [H 

PAIRITEL: Near-infrared observations in the J, .ff and K s were obtained follow ing the VHE 



flare w ith the 1.3m Peters Automated Infrared Imaging Telescope (PAIRITEL; see iBloom et al 



(|2006l )) located at the Fred Lawrence Whipple Observatory. The resulting light curves using dif- 



ferential photometry with four nearby calibration stars are shown in Figured! 



Keck: The optical spectrum of 3C 66A was measured with the LRIS spectrometer (|Oke et al 



19951 ) on the Keck I telescope on the night of 2009 September 17 under good conditions. The in- 
strument configuration resulted in a full- width-half- maximum of ~ 250 km s _1 over the wavelength 
range 3200-5500 A (blue side) and ~ 200 km s" 1 over the ranee 6350-9000A (red side). A series of 
exposures totaling 110 seconds (blue) and 50 seconds (red) were obtained, yielding a signal-to-noise 
(S/N) per resolution element of ~ 250 and 230 for the blue and red cameras respectively. The 
data were reduced with the LowReduj|§ pipeline and calibrated using a spectrophotometric star 
observed on the same night. Inspection of the 3C 66 A spectrum reveals no spectral features aside 
from those imposed by Earth's atmosphere and the Milky Way (Ca H+K). Therefore, these new 
data do not offer any insight on the red shift of 3C 66A an d in particular are unable to confirm the 



previously reported value of z = 0.444 ( Miller et al. 19781 ) 



2.6. Radio Observations 



Radio observations are available thanks to the F-GAMMA (Fermi-Gamma-ray Space Telescope 
AGN Multi-frequency Monito ring Alliance) program, which is dedicated to m onthly monitoring of 



selected Fermi-L AT blazars (Fuhrmann et al. 



2007 



Angelakis et al.l 12008) . Radio flux density 
measurements were conducted with the 100-m Effelsberg radio telescope at 4.85, 8.35, 10.45, and 
14.60 GHz on 2008 October 16. These data are supplemented with an additional measurement 
at 86 GHz conducted with the IRAM 30-m telescope (Pico Veleta, Spain), on 200 8 October 8. 
The data were reduced using standard procedures described in lFuhrmann et all (j2008l ). Additional 
radio observations taken between October 5 and October 15, 2008 (contemporaneous to the flare 
period) are provided by the Medicina, Metsahovi, Noto, and UMRAO observatories, all of which 
are members of the GASP-WEBT consortium. 



ht t p : / / www. ucolick.org/~xavier/ LowRedux /index . html 
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3. Discussion 
3.1. Light Curves 

The resulting multiwavelength light curves from this campaign are shown in Figure [3] for those 
bands with long-term coverage and in Figure U] for those observations that were obtained shortly 
before and after the gamma-ray flare. The VERITAS observations are combined to obtain nightly 
(E > 200 GeV) flux values since no evidence for intra-night variability is observed. The highest 
flux occurred on MJD 54749 and significant variability is observed during the whole interval (x 2 
probability less than 10 -4 for a fit of a constant flux). 

The temporal dependence of the Fermi-LAT photon index and integral flux above 100 MeV 
and 1 GeV are shown with time bins with width of 3 days in Figure O For those time intervals 
with no significant detection a 95%-confidence flux upper limit is calculated. The flux and photon 



index from the Fermi-LAT first source catalog (jAbdo et al.ll2010al ) are shown as horizontal lines for 



comparison. These values correspond to the average flux and photon index measured during the 
first eleven months of Fermi operations, and thus span the time interval considered in the figures. 
It is evident from the plot that the VHE flare detected by VERITAS starting on MJD 54740 is 
coincident with a period of high flux in the Fermi energy band. The photon index during this time 
interval is consistent within errors with t he average photon index T = 1.95 ± 0.03 measured during 



the first 6 months of the Fermi mission (jAbdo et al.ll2010bl ). 



Long-term and well-sampled light curves are available at optical and near-infrared wavelengths 
thanks to observations by GASP-WEBT, ATOM, MDM and PAIRITEL. Unfortunately radio 
observations were too limited to obtain a light curve and no statement about variability in this 
band can be made. The best sampling is available for the R band, for which variations with a 
factor of > 2 are observed in the long-term light curve. Furthermore, variability o n time scales of 



l ess t han a day is observed, as indicated in Figure HI and as previously reported bv lBottcher et al 



(2009) following the WEBT (Whole Earth Blazar Telescope) campaign on 3C 66A in 2007-2008. 



The increase in gamma-ray flux observed in the Fermi band seems contemporaneous with a 
period of increa sed flux in the optical, a nd to test this hypothesis, the discrete correlation function 



(DCF) is used (jEdelson fc Krolik Ill988h . Figure [5] shows the DCF of the F(E > 1 GeV) gamma- 
ray band with respect to the R band with time-lag bins of 3, 5 and 7 days. The profile of the 
DCF is consistent for all time-lag bins, indicating that the result is independent of bin size. The 
DCF with time-lag bins of 3 days was fitted with a Gaussian function of the form DCF(t) = 
C m ax x exp (r — To) 2 /cr 2 , where C max is the peak value of the DCF, to is the delay timescale at 
which the DCF peaks, and a parametrizes the Gaussian width of the DCF. The best fit function 
is plotted in Figure [5] and the best fit parameters are C max = 1.1 ± 0.3, tq = (0.7 ± 0.7) days and 
a = (3.3 ± 0.7) days. An identical analysis was also performed between the F(E > 100 MeV) and 
the R optical band with consistent results. This indicates a clear correlation between the Fermi- 
LAT and optical energy bands with a time lag that is consistent with zero and not greater than ~5 
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Fig. 3. — 3C 66A light curves covering 2008 Aug 22 to 2008 Dec 31 in order of increasing wavelength. The 
VERITAS observations are combined to obtain nightly flux values and the dashed and dotted lines represent the 
average flux measured from the 2007-2008 data and its standard deviation. The Fermi-LAT light curves contain time 
bins with a width of 3 days. The average flux and average photon index measured by Fermi-LAT during the first 
six months of science operations are shown as horizontal lines in the respective panels. In all cases the Fermi-LAT 
photon index is calculated over the 100 MeV to 200 GeV energy range. The long-term light curves at optical and 
infrared wavelengths are presented in the two bottom panels. In the bottom panel GASP-WEBT and PAIRITEL 
observations are represented by open and solid symbols, respectively. 



-19- 



c-" °- 6 f 

\> 0.5 

'e 0.4 ; 

o 

0.3 I 
° 0.21 
",.0.1 

^ 8.0 

E 6.0 

o 

5 2.0 

L_ 

p 8.o : 
« 

E 6.0 : 

o 

a> 

4.0 ; 
° 2.0 : 



VERITAS 

(E>200 GeV) 



Fermi-LAT 

(E > 1 GeV) 



Swift - XRT 

(2-10 keV) 



10 ; 

8 ; 

5 el 
1 4l 



— i 1 1 1 — 

Swift-UVOT/MDM 

V(+1mJy) 
B 

U(-1mJy) 

H 1 1 1 



i 



* i 



Swift - UVOT 

• UVW1 (+2 mJy) 
■ UVM2 

▲ UVW2 (-2 mJy) 



[] 



a 



A A *4 # 



I 1 I 
I I 



' 1 1- 

OGASP - R 
A ATOM - R 
□ MDM - R 



Optical g 



CD SB 03 



-+- 



20 



PAIRITEL 

• H (+5 mJy) 
■ J 

A Ks (-5 mJy) 



740 



745 

MJD - 54000 [days] 



750 



Fig. 4. — 3C 66A light curves covering the period centered on the gamma-ray flare (2008 Oct 1 - 10). The 
VERITAS and Fermi-LAT panels were already described in the caption of Figure [3] Swift Target-of-Opportunity 
(ToO) observations (panels 3-5 from the top) were obtained following the discovery of VHE emission by VERITAS 
( Swordv 2008h . Swift-UVOT and MDM observations are represented by open and solid symbols, respectively. The 
optical light curve in panel 6 from the top displays intra-night variability. An example is identified in the plot, when 
a rapid decline of the optical flux by AF/At ~ —0.2 mJy hr _1 is observed on MJD 54747. 
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Fig. 5. — Discrete correlation function (DCF) of the F(E > 1 GeV) gamma-ray light curve with respect to the R 
band light curve. A positive time lag indicates that the gamma-ray band leads the optical band. Different symbols 
correspond to different bin sizes of time lag as indicated in the legend. The profile of the DCF is independent of bin 
size and is well described by a Gaussian function of the form DCF(t) = C max x exp (r — To ) 2 /ct 2 . The fit to the 
3-day bin size distribution is shown in the plot as solid black line and the best-fit parameters cire Crnax — 1-1 i 0.3, 
r = (0.7 ± 0.7) days and a = (3.3 ± 0.7) days. 
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days. Despite the sparsity of the VERITAS light curve (due in part to the time periods when the 
source was not observable due to the full Moon) the DCF analysis was also performed to search for 
correlations with either the Fermi-LAT or optical data. Apart from the overall increase in flux, no 
significant correlations can be established. The onset of the E > 200 GeV flare seems delayed by 
about ~5 days with respect to the optical-GeV flare but given the coverage gaps no firm conclusion 
can be drawn (e.g., the flare could have been already underway when the observations took place). 
No such lag is expected from the homogeneous model described in the next section but could arise 
in models with complex energy stratification and geometry in the emitting region. 



3.2. SED and Modeling 

The br oadband SED derived from these observations is presented in Figure[6]and modeled using 



the code of lBottcher Chiang! (|2002l ). In this model, a power-law distribution of ultrarelativistic 
electrons and/or pairs with lower and upper energy cutoffs at 7 m ; n and 7 max > respectively, and 
power-law index q is injected into a spherical region of co-moving radius Rb- The injection rate is 
normalized to an injection luminosity L e , which is a free input parameter of the model. The model 
assumes a temporary equilibrium between particle injection, radiative cooling due to synchrotron 
and Compton losses, and particle escape on a time t esc = Vesc Rb/c, where ?7esc is a scale parameter 
in the range ~ 250-500. Both the internal synchrotron photon field (SSC) and external photon 
sources (EC) are considered as targets for Compton scattering. The emission region is moving with 
a bulk Lorentz factor T along the jet. To reduce the number of free parameters, we assume that 
the jet is oriented with respect to the line of sight at the superluminal angle so that the Doppler 
factor is equal to D = (T [1 — /3 cos6> b s ]) -1 = T, where o b s is the angle of the jet with respect 
to the line of sight. Given the uncertainty in the redshift determination of 3C 66A, a range of 
plausible redshifts, namely z = 0.1,0.2,0.3, and the generally used catalog value z = 0.444 are 
consi dered for the mode l ing. A ll model fits include EBL absorption using the optical depth values 



from iFranceschini et al.1 ()2008l ). 



Most VHE blazars known to date are high synchrotron peaked blazars (HSPs), whose SEDs can 
often be fit satisfactorily with pure SSC models. Since the transition from HSP to ISP is continuous, 
a pure SSC model was fit first to the radio through VHE gamma-ray SED. Independently of the 
model under consideration, the low-frequency part of the SED (< 10 20 Hz) is well fit with a 
synchrotron component, as shown in Figure [H For clarity, only the high-frequency range is shown 
in Figures [7] and El where the different models are compared. As can be seen from the figure, a 
reasonable agreement with the overall SED can be achieved for any redshift in the explored range. 
The weighted sum of squared residuals has been calculated for the Fermi-LAT and VERITAS flare 
data (8 data points in total) in order to quantify the scatter of the points with respect to the model 
and is shown in Table HI The best agreement is achieved when the source is located at z ~ 0.2 — 0.3. 
For lower redshifts, the model spectrum is systematically too hard, while at z = 0.444, the model 
spectrum is invariably too soft as a result of EBL absorption. It should be noted that the EBL 
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Fig. 6. — Broadband SED of 3C 66A during the October 2008 multiwavelength campaign. The observation that 
corresponds to each set of data points is indicated in the legend. As an example, the EBL-absorbed EC+SSC model 
for z = 0.3 is plotted here for reference. A description of the model is provided in the text. 
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model of lFranceschini et al.l (120081) predicts some of the lowest optical dep th values in comparison to 
other models (|Finke et al.ll2010i ; iGilmore et al.ll2009l ; IStecker et al.ll2006l ). Thus, a model spectrum 
with redshift of 0.3 or above would be even harder to reconcile with the observations when using 
other EBL models. 

A major problem of the SSC models with z > 0.1 is that Rb is of the order of > 5xl0 16 cm. 
This does not allow for variability time scales shorter than <1 day, which seems to be in contrast 
with the optical variability observed on shorter time scales. A smaller Rb would require an increase 
in the electron energy density (with no change in the magnetic field in order to preserve the flux level 
of the synchrotron peak) and would lead to internal gamma-gamma absorption. This problem could 
be mitigated by choosing extremely high Doppler factors, D > 100. However, thes e are significantly 



large r than the values inferred from VLBI observations of Fermi-LAT blazars (jSavolainen et al 
201ol M Moreover, all SSC models require very low magnetic fields, far below the value expected 
from equipartition (e# = Lb/L c ~ 10~ 3 << 1), where Lb is the Poynting flux derived from the 
magnetic energy density and L e is the energy flux of the electrons propagating along the jet). Table 
H] lists the parameters used for the SSC models displayed in Figure El 

Subsequently, an external infrared radiation field with ad-hoc properties was included as a 
source of photons to be Compton scattered. For all SSC + EC models shown in Figure (H the 
peak frequency of the external radiation field is set to u ext = 1.4 x 10 14 Hz, corresponding to near- 
IR. This adopted value is high enough to produce E > 100 GeV photons from inverse Compton 
scattering off the synchrotron electrons and at the same time is below the energy regime in which 
Klein-Nishina effects take place. Although the weighted sum of squared residuals for EC+SSC 
models are generally worse than for pure SSC models, reasonable agreement with the overall SED 
can still be achieved for redshifts z < 0.3. Furthermore, all SSC + EC models are consistent 
with a variability time scale of Ai var ~ 4 hr. This is in better agreement with the observed 
variability at optical wavelengths than the pure SSC interpretation. Also, while the SSC + EC 
interpretation still requires sub-equipartition magnetic fields, the magnetic fields are significantly 
closer to equipartition than in the pure SSC case, with Ls/L e ~ 0.1. The parameters of the SSC 
+ EC models are listed in Table [5j 

Models with and without EC component yield the best agreement with the SED if the source is 
located at a redshift z ~ 0.2 - 0.3. Of course, this depe nds on the EBL model use d in the analysis. 
An EBL model that predicts higher attenuation than iFranceschini et al.l (|2008l ) would lead to a 
lower redshift range and make it even more difficult to have agreement between the SED models 
and the data when the source is located at redshifts z > 0.4. Finally, it i s worth men t ionin g that 
the redshift range z ~ 0.2 - 0.3 is in agreement with previous estimates by lFinke et al.1 (120081 ). who 
estimate the redshift of 3C 66A to be z = 0.321 based on the magnitude of the host galaxy, and by 



9 As a caveat, jet models with a decelerating flow JGeorganopoulos fc Kazanasl 20031 : Piner et al. 2008h or with 
inhomogeneous transverse structure (Ghisel lini et al . 2005; Irlenri fc Saugefeood l can accommodate very high Doppler 
factors in the gamma-ray emitting region and still be consistent with the VLBI observations of the large scale jet. 
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Fig. 7. — SSC models for redshifts 2 =0.444, 0.3, 0.2, and 0.1 from top to bottom. The Fermi-LAT and VERITAS 
data points follow the same convention used in Figures [T] and [6] to distinguish between flare (red) and dark run (blue) 
data points. In each panel the EBL-absorbed model is shown as a solid red line and the de-absorbed model as a red 
dashed line. De-absorbed VERITAS flare points are shown as open squares. In all cases the optical depth values from 
Franceschini et al, ( 20081) are used. The best agreement between the model and the data is achieved when the source 



is located at z = 0.2 — 0.3. For lower redshifts the model spectrum is systematically too hard, while at z — 0.444, the 
model spectrum is too soft. 
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Fig. 8. — EC + SSC model for redshifts z =0.444, 0.3, 0.2, and 0.1 from top to bottom. The individual EBL- 
absorbed EC and SSC components are indicated as a dash-dotted and dotted lines, respectively. The sum is shown 
as solid red line (dashed when de-absorbed). The best agreement between the model and the data is achieved when 
the source is located at z ~ 0.2. 
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Prandini et al.l (|201Cl ) who use an empirical relation between the previously reported Fermi-LAT 
and IACTs spectral slopes of blazars and their redshifts to estimate the redshift of 3C 66A to be 
below z = 0.34 ± 0.05. 

A detailed study of hadronic versus leptonic modeling of the October 2008 data will be pub- 
lished elsewhere, but it is worth mentioning that the synchrotron proton blazar (SPB) model has 
been used to adequately re produce the quasi-s imultaneous SED observed during the 2003-2004 
multi- wavelength campaign (jReimer et al.ll2008l ). On that occasion rapid intraday variations down 
to 2 hours time scale were observed, while during the 2008 campaign presented here these varia- 
tions seem less rapid. Qualitatively, the longer time scale variations may be due to a lower Doppler 
beaming, at the same time that a strongly reprocessed proton synchrotron component dominates 
the high energy output of this source. 



4. Summary 

Multiwavelength observations of 3C 66A were carried out prompted by the gamma-ray outburst 
detected by the VERITAS and Fermi observatories in October 2008. This marks the first occasion 
that a gamma-ray flare is detected by GeV and TeV instruments in comparable time scales. The 
light curves obtained show strong variability at every observed wavelength and in particular, the 
flux increase observed by VERITAS and Fermi is coincident with an optical outburst. The clear 
correlation between the Fermi-LAT and R optical light cur yes permits one to g o beyond the source 



association reported in the 1st Fermi-LAT source Catalog (jAbdo et al.ll2010al ) and finally identify 
the gamma-ray source 1FGL J0222. 6+4302 as blazar 3C 66A. 

For the modeling of the overall SED a reasonable agreement can be achieved using both a pure 
SSC model and an SSC + EC model with an external near-infrared radiation field as an additional 
source for Compton scattering. However, the pure SSC model requires (a) a large emission region, 
which is inconsistent with the observed intra-night scale variability at optical wavelengths and 
(b) low magnetic fields, about a factor ~ 10~ 3 below equipartition. In contrast, an SSC + EC 
interpretation allows for variability on time scales of a few hours, and for magnetic fields within 
about an order of magnitude of, though still below, equipartition. It is worth noting that the results 
presented here agree wi th the findings follow ing the (E > 200 GeV) flare of blazar W Comae (also 



an ISP) in 2008 March (lAcciari et al.ll2008al ). In both cases the high optical luminosity is expected 



to play a key role in providing the seed population for IC scattering. 

Intermediate synchrotron peaked blazars like 3C 66A are well suited for simultaneous obser- 
vations by Fermi-LAT and ground-based IACTs like VERITAS. Relative to the sensitivities of 
these instruments, ISPs are bright enough to allow for time-resolved spectral measurements in each 
band during flaring episodes. These types of observations coupled with extensive multi-wavelength 
coverage at lower energies will continue to provide key tests of blazar emission models. 
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Model parameter 


z = 0.1 


z = 0.2 


z = 0.3 


z = 0.444 


Low-energy cutoff (7min) 


1.8xl0 4 


2.0xl0 4 


2.2xl0 4 


2.5xl0 4 


High-energy cutoff (7 max ) 


3.0xl0 5 


4.0xl0 5 


4.0xl0 5 


5.0xl0 5 


Injection index (q) 


2.9 


2.9 


3.0 


3.0 


Injection luminosity (L e ) [10 45 erg s -1 ] 


1.3 


3.3 


5.7 


12.8 


Co-moving magnetic field (is) [GJ 


0.03 


0.02 


0.02 


0.01 


Poynting flux (L B ) [10 42 erg s -1 ] 


1.1 


4.9 


8.5 


13.7 


es = Lb/L e 


0.9xl0" 3 


1.5xl0~ 3 


1.5xl0~ 3 


l.lxl0~ 3 


Doppler factor (D) 


30 


30 


40 


50 


Plasmoid radius (Rb) [10 16 cm] 


2.2 


6.0 


7.0 


11 


Variability time scale (Ji™ 1 ") [hr] 


7.4 


22.1 


21.1 


29.4 


Weighted sum of squared residuals to VERITAS flare data 


7.1 


0.9 


0.7 


6.2 


Weighted sum of squared residuals to Fermi-LAT flare data 


1.6 


1.6 


1.3 


1.4 


Total weighted sum of squared residuals 


8.7 


2.5 


1.9 


7.6 



Table 4: Parameters used for the SSC models displayed in Figure[71 All SSC models require very low 
magnetic fields, far below the value expected from equipartition (i.e. e# << 1). The weighted sum 
of squared residuals to the VERITAS and Fermi-LAT data and the total value for the combined 
data set are included at the bottom of the table. The best agreement between the model and the 
data is obtained when the source is at redshift z = 0.2 — 0.3. 



Model parameter 


z = 0.1 


z = 0.2 


z = 0.3 


z = 0.444 


Low-energy cutoff (7 m in) 


5.5xl0 3 


7.0xl0 3 


6.5xl0 3 


6.0 xlO 3 


High-energy cutoff (7 max ) 


1.2 xlO 5 


1.51.2xl0 5 


1.51. 2xl0 5 


1.51. 2xl0 5 


Injection index (q) 


3.0 


3.0 


3.0 


3.0 


Injection luminosity (L e ) [10 44 erg s~ 4 ] 


1.1 


4.2 


6.0 


10.4 


Co-moving magnetic field (B) [G] 


0.35 


0.22 


0.21 


0.23 


Poynting flux (L B ) [10 43 erg s _1 ] 


1.0 


2.4 


6.0 


11.2 


6b = Ls/L e 


0.10 


0.06 


0.10 


0.11 


Doppler factor {D) 


30 


30 


40 


50 


Plasmoid radius {Rb) [10 16 cm] 


0.5 


1.2 


1.5 


1.5 


Variability time scale (8t™™) [hr] 


1.7 


4.4 


4.5 


4.0 


Ext. radiation energy density [10 -6 erg cm~ 3 ] 


5.4 


2.4 


1.2 


1.3 


Weighted sum of squared residuals to VERITAS flare data 


4.8 


3.6 


7.9 


15.7 


Weighted sum of squared residuals to Fermi-LAT flare data 


1.0 


1.2 


0.8 


1.5 


Total weighted sum of squared residuals 


5.8 


4.8 


8.7 


17.2 



Table 5: Parameters used for the EC+SSC model fits displayed in Figure [HJ These model fits 
require magnetic fields closer to equipartition and allow for the intra-night variability observed in 
the optical data. The weighted sum of squared residuals to the VERITAS and Fermi-LAT data 
and the total value for the combined data set are included at the bottom of the table. 
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